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ABSTRACT 

We present NuSTAR high-energy X-ray observations of the pulsar wind nebula (PWN)/supernova remnant 
G2 1.5— 0.9. We detect integrated emission from the nebula up to ~40 keV, and resolve individual spatial features 
over a broad X-ray band for the first time. The morphology seen by NuSTAR agrees well with that seen by 
XMM-Newton and Chandra below 10 keV. At high energies, NuSTAR clearly detects non-thermal emission up to 
~20 keV that extends along the eastern and northern rim of the supernova shell. The broadband images clearly 
demonstrate that X-ray emission from the North Spur and Eastern Limb results predominantly from non-thermal 
processes. We detect a break in the spatially integrated X-ray spectrum at ~9 keV that cannot be reproduced by 
current spectral energy distribution models, implying either a more complex electron injection spectrum or an 
additional process such as diffusion compared to what has been considered in previous work. We use spatially 
resolved maps to derive an energy-dependent cooling length scale, L(E) oc E m with m = —0.21 zb 0.01. We find 
this to be inconsistent with the model for the morphological evolution with energy described by Kennel & Coroniti. 
This value, along with the observed steepening in power-law index between radio and X-ray, can be quantitatively 
explained as an energy-loss spectral break in the simple scaling model of Reynolds, assuming particle advection 
dominates over diffusion. This interpretation requires a substantial departure from spherical magnetohydrodynamic, 
magnetic-flux-conserving outflow, most plausibly in the form of turbulent magnetic-field amplification. 

Key words: ISM: individual objects (G2 1.5— 0.9) - ISM: supernova remnants - radiation mechanisms: general - 
stars: neutron - X-rays: ISM 
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1. INTRODUCTION 

A pulsar wind nebula (PWN) is a bubble of relativistic 
particles and magnetic field inflated by a rotation-powered 
pulsar, emitting centrally peaked synchrotron radiation. Young 
PWNe are frequently found inside shell supernova remnants 

(SNRs), where the relativistic wind of electron-positron pairs 

(and perhaps ions) experiences a wind termination shock close 
to the pulsar due to the pressure of the SNR interior. Beyond 
this wind shock, the relativistic fluid can radiate synchrotron 

emission from the radio to the X-ray band, and inverse-Compton 

(IC) emission at higher energies. A PWN/SNR combination is 

often called a “composite” SNR. PWNe can outlive their SNR 

and interact directly with the interstellar medium. See Gaensler 

& Slane (2006) for a review. 

The detailed characterization of PWN spectra gives informa- 
tion on the particle energy distribution produced by the pulsar, 

and also on the nature of acceleration in relativistic shocks. How- 
ever, spectral structure present immediately downstream of the 
wind shock can be altered by propagation effects including dif- 
fusive transport and radiative losses, which may depend on the 


evolution of the entire PWN. The ability to use PWNe as labora- 
tories in which to study the behavior of relativistic pair plasmas 
depends largely on the extent to which these various effects can 
be disentangled. Models studying PWNe with time-dependent, 
one-zone, homogeneous approximations, provide insights into 
the evolution of the spectrum across all energy bands of the neb- 
ula as a whole. Broadband spectral energy distributions (SEDs), 
coupled with spatially resolved spectroscopy, are required for 
this purpose. 

Radially dependant models also provide valuable insights by 
addressing detailed spatial and spectral structures of PWNe. The 
classic work of Kennel & Coroniti (1984), which invokes a par- 
ticular hydrodynamic model, predicts the spectral break between 
the optical and X-ray energy bands, as well as the behavior of 
size with photon energy for the Crab Nebula. Implicit in that 
calculation is a prediction for the radial dependence of the spec- 
trum. Kennel & Coroniti (1984) assume particle transport by 
pure advection in a spherical geometry; this situation predicts a 
roughly uniform nebular spectrum with radius until the (energy- 
dependent) nebular edge, where radiative energy losses sharply 
steepen the spectrum. This behavior is generally inconsistent 
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with observations (Reynolds 2003; Tang & Chevalier 2012), in 
particular in several objects a fairly uniform spectral steepening 
with radius is observed (see Bocchino & Bykov 2001). 

The cause of spectral structure is of interest beyond the spe- 
cific PWN context; if its origin is in the physics of particle 
acceleration at relativistic shocks, the results have implications 
for other objects which share the same shock mechanisms. How- 
ever, the structure may be due instead to transport and evolu- 
tionary effects in the PWN post-shock flow. At the same time, 
spatially integrated spectra of PWNe are not well understood, 
especially in their most salient feature — a large steepening of 
the spectrum between the flat radio emission and the consider- 
ably steeper spectra observed from IR through X-rays. Chevalier 
(2005) documents differences between radio and X-ray spectral 
indices of A = a x — a r ~ 0.7-1 (S v a v~ a ) for seven out of 
the eight PWNe modeled, compared to radiative energy losses 
in homogeneous steady sources which can produce, at most, 
A = 0.5. If this break is not due solely to evolutionary effects 
and instead is due to at least, in part, to intrinsic spectral struc- 
ture in the particle distribution injected at the wind shock, it 
is implying something important and interesting about particle 
acceleration in pulsar winds. 

G21.5— 0.9 was discovered in 1970 (Altenhoff et al. 1970; 
Wilson & Altenhoff 1970) in the radio band, and first observed 
in the X-ray band in 1981 (Becker & Szymkowiak 1981). 
It is a classic example of a Crab-like PWN: it has a filled, 
mostly symmetric spherical morphology centered on a pulsar. 
Measurements spanning 44 yr show a well-characterized flat 
spectrum in the radio regime (e.g., Goss & Day 1970; Becker 
& Kundu 1976; Salter et al. 1989; Bandiera et al. 2001), with 
recent observations reporting a spectral index of a r = 0.0 zb 0. 1 
(Bietenholz et al. 2011). The flux density at 1 GHz is 6 Jy 
(Camilo et al. 2006), and 13 CO and Hi analyses determined the 
distance to the nebula to be 4.7 =b 0.4 kpc (Camilo et al. 2006; 
Tian & Leahy 2008). In this paper we adopt a distance of 5 kpc. 

Observations with Chandra and XMM-Newton find that the 
X-ray emission is dominated by a centrally peaked core that 
contains ~85% of the 2-8 keV flux. The spectrum is described 
by a non-thermal power law with no evidence of line emission. 
The total unabsorbed nebular flux is F x (0.5-10 keV) = 9.35 x 
10 -11 erg cm -2 s -1 (Slane et al. 2000; Safi-Harb et al. 2001; 
Warwick et al. 2001). Safi-Harb and coworkers detected spectral 
steepening in the nebula, indicative of synchrotron cooling. The 
innermost O'/ 5 radius region has a power-law photon index of 
r = 1.43 ± 0.02 (a x = T— 1) which softens to T = 2.13 ±0.06 
at a radius of 40" consistent with the edge of the nebula. 
This spectral softening is also visible in hardness ratio images 
(Matheson & Safi-Harb 2005). The photon index of T ~ 2 
yields A ~ 1 . 

The associated pulsar, PSR J1833— 1034, was discovered 
in 2005, 35 yr after the initial detection of the PWN. PSR 
J1833— 1034 has a 61.86 ms period, P = 2.0 x 10 -13 , 
r c = 4.8 kyr, and E = 3.3 x 10 37 erg s -1 (Gupta et al. 
2005; Camilo et al. 2006). Despite the pulsar’s high spin- 
down luminosity, pulsations have not been detected in the 
X-ray band. Chandra and XMM-Newton observations detected 
a region of diffuse, uniform emission extending from the edge 
of the PWN at 40" out to 150" (e.g., Slane et al. 2000), with 
a non-thermal power-law spectrum of T ~ 2.5. The flux is 
substantially dimmer than the PWN. This symmetric emission 
was proposed to be an extension of the PWN itself (Warwick 
et al. 2001), but the absence of coincident radio emission, and 
the recognition of the importance of dust scattering for the large 


column density toward G2 1.5— 0.9, led Bandiera & Bocchino 
(2004) and Bocchino (2005) to model the smoothly distributed 
halo emission as due to dust scattering. 

Additionally, two regions of brightened emission were dis- 
covered near the outer edges of the halo, referred to as the North 
Spur and Eastern Limb (Warwick et al. 2001; Safi-Harb et al. 
2001) that could not be explained by dust scatter. These re- 
gions, located approximately 80" and 120", respectively, from 
the center of the nebula, have weak but detectable X-ray emis- 
sion (Bocchino 2005; Bocchino et al. 2005; Matheson & Safi- 
Harb 2010). Deep radio observations have not detected any 
emission from either the halo or the Eastern Limb, while the 
North Spur can be clearly seen in a 1.4 GHz image (Bietenholz 
et al. 2011). Bocchino (2005), and later Matheson & Safi-Harb 
(2010), found the spectrum of the North Spur comprises of a 
weak, low-temperature thermal component and a non-thermal 
continuum. This knot of emission has been interpreted as the 
result of ejecta interacting with the H-envelope of the supernova 
(SN). The Eastern Limb, in contrast, was found to have only 
a non-thermal spectrum, and has spectral and morphological 
features that imply it is the limb-brightened region of the SN 
shell (Matheson & Safi-Harb 2010). However, the low surface 
brightness of both the Eastern Limb and North Spur prevent their 
<10 keV continua from being characterized more specifically. 

In this paper we present the first subarcminute X-ray images 
above 10 keV and corresponding X-ray spectroscopic studies 
of G21.5— 0.9. Section 2 discusses the NuSTAR observations. 
Section 3 presents our spectral analysis, while Section 4 presents 
our image analysis and Section 5 details our search for the pulsar 
in the high-energy X-ray band. Lastly, in Section 6 we discuss 
the spectral and spatial studies and how NuSTAR can shed light 
on the physics in a PWN and the natures of the North Spur and 
Eastern Limb. 

2. NuSTAR OBSERVATIONS 

NuSTAR observed G21.5— 0.9 on four separate occasions 
for a total of 281 ks: 2012 July 29 (ObsID 10002014003), 
2012 July 30 (ObsID 10002014004), 2013 February 26 (ObsID 
40001016002), and 2013 February 27 (ObsID 40001016003). 
The center of the remnant was located approximately 2' from 
the on-axis position for all observations so that the majority 
of the PWN was located on one of the four detector chips. 
NuSTAR , which contains two co-aligned optic/detector focal 
plane modules (FPMA and FPMB), has a half-power diameter 
of 58", an angular resolution of 18" (FWHM) over its 3-79 keV 
X-ray energy range, and a characteristic FWHM energy res- 
olution of 400 eV at 10 keV. The field of view is 12' x 12' 
at 10 keV as defined by the full width at half-intensity. 
The NuSTAR nominal reconstructed coordinates are accurate 
to 8" (90% confidence level; Harrison et al. 2013). 

Prior to our imaging and spectral analysis, we registered 
individual observations to J2000 coordinates using the central 
peak position of the PWN measured by Chandra : R.A.(J2000) = 
18 h 33 m 33 s .54, decl.(J2000) = -10°34'07"6 (Safi-Harb et al. 
2001). Note that the Chandra centroid position is <1" offset 
from the radio pulsar position (Camilo et al. 2006). The 
NuSTAR field of view for G2 1.5— 0.9 is devoid of any visible 
point sources including the foreground star SS 397, which 
is located ~100" southwest of the center of the PWN. We 
determined the centroid position of the inner 30" radius region 
of the PWN in the 3-10 keV band using the IDL routine 
gcntrd so that the NuSTAR image is not contaminated by 
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substructures such as the Eastern Limb and North Spur. The 
centroiding errors (90% confidence level) are ~3" in both R.A. 
and decl. We subsequently confirmed by fitting a circular two- 
dimensional (2D) Gaussian convolved with the NuSTAR point- 
spread function (PSF; Section 4.1) that the centroid position 
of the PWN overlaps with the radio pulsar position within our 
uncertainty. 

3. SPECTROSCOPY 

We performed NuSTAR spectral analysis integrated over the 
PWN region as well as spatially resolved spectroscopy which 
we present in subsequent sections. We applied the same analysis 
procedures, described below, for all NuSTAR spectra. 

Prior to spectral fitting we generated NuSTAR response matrix 
(RMF) and effective area (ARF) files for an extended source 
using nuproducts (NuSTARDAS vl.1.1), then grouped the 
spectra to >20 counts per bin using FTOOLS grppha (HEAsoft 
6.13). We fit the NuSTAR spectra using XSPEC version 12.8.0 
(Arnaud 1996), with the atomic cross sections set to those from 
Verner et al. (1996) and the abundances to Wilms et al. (2000). 
We fit the power-law model pegpwrlw and the broken power- 
law model bknpower. We multiplied the continuum models 
by the interstellar absorption model, Tbabs, with Ah frozen 
to 2.99 x 10 22 cm -2 from Tsujimoto et al. (2011). We fit 
the spectra from each observation and each module jointly by 
linking all parameters except the continuum normalization, to 
take into account small calibration uncertainties between the 
two modules’ flux normalization. 

We generated background spectra using Nulyses , NuSTAR- 
specific software that accounts for detector background and 
cosmic X-ray background (CXB) for a given extraction region 
of the source spectrum. Nulyses creates source-free background 
maps from the NuSTAR blank-sky survey data taken less than 
a month from the G2 1.5— 0.9 observations. This method was 
applied to all subsequent spectral analyses. The conventional 
way of extracting background spectra from a nearby region 
is not applicable to our analysis since the PWN, broadened 
by the NuSTAR PSF, covers most of the detector chip and the 
detector background varies among the different chips. We fit the 
NuSTAR spectra in the 3-45 keV band, above which the detector 
background dominates. 

3.1. NuSTAR Spectroscopy of the Entire PWN 

We extracted NuSTAR data using a 165" radius circular 
region centered at the pulsar PSR J1 833— 1034 position. This 
region represents the largest circle within the same detector 
chip (to ensure that NuSTAR detector background is the same 
and predictable), and it encloses over 93% of the PWN photons. 

We first fit the NuSTAR spectra with a single absorbed power- 
law model in the 3-45 keV energy band. We obtained a best- 
fit power-law index of T = 2.04 ±0.01 and an X-ray flux 
of E 2 -8keV = 5.54 ± 0.03 erg s -1 cm -2 . These results are 
consistent within their uncertainties to the parameters reported 
by Tsujimoto et al. (2011) (r = 2.05 ± 0.04, T^-skeV = 
5.7 ± 0.5 erg s -1 cm -2 ), who analyzed G21.5— 0.9 for the 
purposes of cross-calibrating X-ray instruments operational at 
that time. All uncertainties (90% confidence level) quoted in the 
text and tables include both systematic and statistical errors. A 
more detailed comparison of G21 .5 — 0.9 spectral fitting between 
NuSTAR and other X-ray telescopes as well as systematic error 
estimates will be addressed in a separate calibration paper. 

A single power-law fit is not satisfactory, with reduced x 2 — 
1.33, with clear residuals evident between 5 and 10 keV (see the 
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Figure 1 . NuSTAR spectral fitting of G2 1.5— 0.9 in the energy range 3-45 keV 
from the entire PWN (extraction radius = 1 65"). Eight NuSTAR spectra and their 
residuals are shown. The middle panel depicts the residuals for the absorbed 
power-law fit, while the lower panel depicts the absorbed broken power-law fit. 
(A color version of this figure is available in the online journal.) 


left panel of Figure 1). We then fit a broken power-law model 
(bknpower) to the NuSTAR spectra, yielding reduced x 2 = 
1 .09. Statistical tests using the F-distribution (f t est in XSPEC) 
show the broken power-law model is statistically favored over 
the single power-law model with high significance, so that a 
spectral break is unambiguously required. The best-fit power- 
law indices were Ti = 1.996 ± 0.013 and r 2 = 2.093 ± 0.013, 
with the best- fit break energy of £break = 9.7 ± 1.3 keV. 

In order to investigate any instrumental effects that mimic 
a break around 9 keV, we analyzed NuSTAR observations of 
the known power-law source 3C273. We followed the same 
procedures used for G21.5— 0.9, and found no systematic 
residuals around 9 keV in the NuSTAR data. A spectral break 
in another PWN, the Crab Nebula, has recently been detected 
by NuSTAR at similar energies of 8-11 keV (K. C. Madsen, in 
preparation). 

Dust scattering can also mimic the low-energy spectral 
behavior observed in the NuSTAR spectra. Recent analysis of 
two low-mass X-ray binaries with similar absorption columns 
of~3x 10 22 cm -2 ,GX5— 1 and GX 13+1, present radial profiles 
of models of the dust scattering halos. These profiles indicate 
that 0. 1 % of the source photons at 2.5 keV are enclosed between 
50" < r < 600" (Smith et al. 2002). Including the energy- 
dependant halo intensity profile relationship I oc E~ 2 (Predehl 
& Schmitt 1995; Smith 2008), a radius of r < 160" encloses 
over 99% of the source photons in both 3-5 keV and 5-8 keV. 
Similar analysis reveals that the r < 30" extraction region, 
presented in Section 3.2, has a scattering loss of ~0.2% and 
~0. 1% of source photons in 3-5 keV and 5-8 keV, respectively. 
Loss of photons by dust scattering therefore has a negligible 
effect, and cannot explain the spectral break in the NuSTAR 
spectrum. Additionally, because our scattering estimates are 
based on point sources, we cannot conclusively rule out the 
possibility of a very small effect of dust scattering on the radial 
behavior of the low-energy spectrum presented in Section 3.2. 
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Table 1 

Spectral Fitting of the NuSTAR G21.5— 0.9 Data 


Parameter 
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9.74 ± 1.0 


Let!— 8 keV) 

5.45 ± 0.03 
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F x ( 15-50 keV) 

5.47 ± 0.03 

5.11 ±0.08 






X 2 /dof 

1.33(5117) 

1.09(5112) 

1.32(3564) 

0.99 (3557) 

1.08(2257) 

0.960(2255) 

0.942(1924) 


Notes. Spectral fitting of G2 1.5— 0.9 with various extraction regions. The column header indicates the region from which the data was extracted. Ah 
was frozen to 2.99 x 10 22 cm -2 for all fits. Flux is listed in units of 10 “ 11 erg s -1 cm -2 . The goodness of fit is evaluated by the reduced y } . The errors 
are 90% confidence level. 



Figure 2. NuSTAR G21 .5— 0.9 data, represented by red crosses, overplotted on a 
graph containing data from previous observations, represented by black crosses, 
and a model prediction of the spectrum of G2 1.5— 0.9 shown by the solid black 
line. From Tanaka & Takahara (2011). 

(A color version of this figure is available in the online journal.) 


An X-ray spectral break at ~12 keV was first suggested by 
Tanaka & Takahara (2011). They extrapolated the spectral fits 
from the hard and soft X-ray bands reported by Tsujimoto 
et al. (2011), and noted that a break at ~12 keV was likely. 
NuSTAR not only confirms the break, but provides the first 
measurement of the ~9 keV break in broadband 3-45 keV 
spectra of G2 1.5— 0.9 in a single X-ray telescope. 

Obtaining a well-calibrated X-ray spectrum from G2 1.5— 0.9 
is particularly valuable when creating PWN theoretical models. 
Specifically, Tanaka & Takahara (2011) presented synchrotron 
emission from radio to GeV energies by modeling the electron 
spectrum emanating from the pulsar of G2 1.5— 0.9. They 
incorporated adiabatic and energy losses in a one-zone model. 
Figure 2 shows the NuSTAR best-fit broken power-law data 
overplotted on Figure 3 from Tanaka & Takahara (2011). More 
details are presented in Section 6.1. 

3.2. Spatially Resolved Spectroscopy 

A number of young PWNe, including G21.5— 0.9, exhibit 
spectral softening from the center of the PWN outward. This 
is the classic signature of synchrotron burn-off. Slane et al. 
(2000) and later Safi-Harb et al. (2001) extracted Chandra 
spectra at various annuli of G2 1.5— 0.9 and showed spectral 
softening from T = 1.43 ± 0.02 in the central 5" radius circle 
to r = 2.13 zb 0.06 in the outer annulus at a radius of ~40". 



Radial Distance from Pulsar Position [arcsec] 

Figure 3. NuSTAR spectral photon indices as a function of extraction region 
for an absorbed bknpower fit. The radial positions are measured from the 
radio location of the pulsar, R.A.(J2000) = 18 h 33 m 33®54, Decl.(J2000) = 
— 10 o 34'07"6. The blue points illustrate the photon index T\ below the spectral 
break of ~9 keV. The red points illustrate the photon index r 2 above the break. 
(A color version of this figure is available in the online journal.) 


NuSTAR has the ability to probe the synchrotron burn-off effects 
in the hard X-ray band above 10 keV. 

We extracted NuSTAR spectra from the central 30" region 
as well as from four nested annuli, each 30" in width. After 
experimenting with various annulus widths, we found that the 
annulus width of 30" is the best compromise between the 
NuSTAR angular resolution and the expected spatial variation 
from the Chandra results. 

The broadband 3-45 keV fit results from both the absorbed 
power law and bknpower models are shown in Table 1. The 
power-law index fit with a single power-law model increases 
fromT = 1. 97 ± 0.01 in the 30" inner region to T = 2.25±0.02 
in the outer 90"- 120" radius annulus, confirming the spectral 
softening discovered by Chandra. However, the residuals for 
the inner 30" and 30"-60" regions also show a spectral break 
around 9 keV. We fit a broken power-law model, and the results 
are shown in Table 1 and plotted in Figure 3. We conclude that 
the spectral break at 9 keV is detected with high significance in 
the inner 30" radius circular region and the 30"-60" annulus. 
Interestingly, as is shown in Figure 3, the second (high energy) 
power-law component shows no spatial variation. The spectral 
softening is observed only below 9 keV. 

4. IMAGING ANALYSIS 

In the imaging analysis, we applied common procedures for 
image preparation. Specifically, we generated the following: (1) 
mosaic images from the four observations, (2) exposure maps 


4 


The Astrophysical Journal, 789:72 (12pp), 2014 July 1 


Nynka et al. 


with vignetting effects included, and (3) an effective PSF that 
takes into account the various off-axis positions, orientation 
angles, and exposure times of the observations. After proper 
source registration as described in Section 2, we selected photon 
events in different energy bands using dmcopy (CIAO v.4.4). 
We chose energy bands to ensure sufficient photon statistics 
in each image, as well as to minimize the effects of averaging 
the energy-dependent vignetting function over the energy range. 
We then summed all four NuSTAR raw images in each energy 
band to create mosaic images using XIMAGE (HEAsoft 6.13). 
For each observation, we generated an exposure map with the 
optics vignetting effects using nuexpomap (NuSTARDAS v 
1.1.1) and summed those as well. We then applied exposure- 
corrections to the raw mosaic images in XIMAGE and generated 
energy-dependent flux images. Figure 5 shows the 3-6 keV 
flux image on the left, clearly showing the centrally peaked 
PWN, broadened by the wings of the NuSTAR PSF. In order to 
investigate the radial profile and detect subtle features possibly 
buried in the raw images, we generated an effective PSF used for 
image deconvolution and 2D forward image fitting, as described 
in the subsequent sections. 

We estimated the background level using the nulysis software 
described in Section 3. First, we reproduced the background 
levels in the three detector chips where the contribution of 
G2 1.5— 0.9 is negligible. We found that the background for the 
G2 1.5— 0.9 observation is largely due to the stray-light CXB 
component below 20 keV and the internal detector background 
above 20 keV. We were able to reliably produce the background 
count rates on the detector chip containing the G2 1.5— 0.9 
image, and found that the source emission is dominant up to 
30 keV. Hereafter we present imaging analysis below 30 keV. 

4.1. G2 1.5— 0.9 PWN Energy -dependent Radius 

A measurement of the PWN radius in different energy bands 
is complementary to the radially dependent spectral analysis. 
The angular resolution of NuSTAR is comparable to the PWN 
size (~80"). Therefore, it is essential to take into account the 
blurring effects by the PSF convolution. 

We adopted a forward-folding method to measure the size of 
the PWN using Sherpa, CIAO’s modeling and fitting package 
(Fruscione et al. 2006). The use of Sherpa allows us to use the 
goodness-of-fit test (C-statistics) to find the best-fit parameters. 
We convolved an assumed source model, a circular 2D Gaussian 
profile, with the effective PSF described in Section 3, and fit the 
NuSTAR flux images. We also ran the conf command in Sherpa 
to determine 90% confidence intervals for all fit parameters. 

We confined our fitting range to within a radius of 160". 
In general, a circular 2D Gaussian provides a good fit to the 
entire PWN, leaving only small residuals. The best-fit Gaussian 
FWHMs in five energy bands are plotted in Figure 4. The trend 
of decreasing PWN size with energy is evident, confirming 
the synchrotron bum-off effect at energies above 10 keV. In 
Section 6.2 we employ the FWHM radius to measure the 
synchrotron cooling length. We fit the data with a power-law 
L(E) oc E m . The best fit of L(E ) is shown in Figure 4 and 
yields an index of m = — 0.21±0.01. These results are used in 
Section 6.2 to constrain physical conditions in the nebula. 

4.2. The North Spur and Eastern Limb 

The Eastern Limb and North Spur are only partially resolved 
in the raw NuSTAR mosaics due to the small size of G21.5— 0.9 
compared to the NuSTAR PSF, combined with the low surface 


FWHM as a function of Energy 



Figure 4. 2D Gaussian FWHM radius of G21.5— 0.9 as a function of energy. 
The midpoint in each band is the mean energy, weighted by the NuSTAR flux. 
The data is well fit to a power-law model L(E) oc E m , with m = —0.21 ± 0.01. 
(A color version of this figure is available in the online journal.) 

brightness of these features. To detect these features with 
high confidence, we applied two different methods. First, we 
confirmed the detection of the Eastern Limb and Northern 
Spur in the flux images in each energy band (Section 4.2.1). 
Second, we sharpened the NuSTAR flux images using the 
Lucy-Richardson deconvolution algorithm (Richardson 1972; 
Lucy 1974). 

4.2.1. ID Profile Analysis 

In order to detect the Eastern Limb and North Spur in the 
raw NuSTAR image, we analyzed line intensity profiles of 
6 image pixels (a total of ~9") in width. We chose projection 
axis lines along different orientation angles around the center 
of the PWN. The left image in Figure 5 contains two lines, 
colored red and blue, that indicate the positions along which the 
flux profiles were taken. The red line was chosen to fall along 
the regions of brightest Eastern Limb emission, while the blue 
line overlays a region of G2 1.5— 0.9 that does not contain any 
emission from either the Eastern Limb or the North Spur. The 
intensity profiles on the right of Figure 5 correspond to the lines 
drawn in the left of Figure 5. The red profile shows clear and 
significant excess emission when compared to the blue profile, 
confirming that the Eastern Limb is detected in the 3-6 keV 
band. This also holds for the 6-10 keV and 10-15 keV. While 
the two profiles in the 15-20 keV band are very similar, the lines 
are separate and distinct between 60"-40", when accounting for 
statistical uncertainty. With these intensity profiles we confirm 
the detection of the North Spur and Eastern Limb at energies as 
high as 20 keV. 

We repeated this process with the North Spur, orienting the 
red line along the north/ south axis. Lastly, we also confirmed the 
detection of the Eastern Limb and North spur in residual maps 
from fitting the PWN, as mentioned in Section 4. 1 . The residual 
maps had areas of faint excess, indicating that there exists 
emission aside from the PWN. However, the broad NuSTAR 
PSF prevents us from studying any further morphology of the 
faint emission. We therefore turn to image deconvolution in an 
attempt to remove the effects of the PSF. 

4.2.2. Image Deconvolution: Method and Verification 

We applied an iterative deconvolution technique to the NuS- 
TAR images using are store (CIAO v4.4) and the effective 
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Figure 5. Intensity profiles of G2 1.5— 0.9 were obtained to confirm the existence of the Eastern Limb and North Spur in the raw NuSTAR images. Left: 3-6 keV 
NuSTAR mosaic image. Exposure-map vignetting corrections were applied, and FPMA and FPMB summed. The red and blue lines indicate the locations along which 
the profiles were obtained. Right: intensity profiles as a function of distance (in arcseconds) from the PWN center from the NuSTAR 3-6 keV, 6-10 keV, 10-15 keV, 
15-20 keV image. The profiles show clear significant excess across the north eastern side of G21.5— 0.9, confirming the detection of excess emission. 

(A color version of this figure is available in the online journal.) 



Figure 6. Deconvolved NuSTAR images at various energy bands: (b) 3-6 keV, (c) 6-10 keV, (d) 10-15 keV, (e) 15-20 keV, and (f) 20-25 keV. The images show the 
faint emission from the Eastern Limb and North Spur. The images are shown on a logarithmic scale, and colors were chosen to highlight the non-plerionic details. 
Image (a) shows the Chandra 3-6 keV image for comparison. The green circle has a radius of 165". 

(A color version of this figure is available in the online journal.) 


PSFs described in Section 4. The iterative image deconvolu- 
tion can produce artificial features if the process is over-iterated 
and/or the background region is deconvolved. Great care was 
taken to ensure that the deconvolved image was both stable and 
reproducible. We deconvolved the 3-6 keV mosaiced NuSTAR 
image with several iterations, from 20 up to 200 in increments 
of 20. Each of these images exhibited the same features of the 
North Spur and Eastern Limb that are characterized by areas 
of brightened emission to the north and northeast of the PWN. 
Additionally, we confirmed that the features visible through de- 
convolution were not dependent on telescope rotation, detector 
module, or observation. The data sets were first grouped by 


NuSTAR module, summing all the FPMA images before decon- 
volving, and likewise for all the FPMB images. The deconvolved 
images have identical features. We also grouped the data by 
epoch, summing the 2012 and 2013 data separately. As before, 
the two deconvolved images are very similar. 

As a final verification of our deconvolution process, we 
compared the output images to the 3-6 keV image from 
Chandra. Figure 6(a) shows the Chandra image at the top left 
corner. One can see the Eastern Limb and North Spur, which is 
identified as “knot” in the image. The North Spur is visible as the 
excess of emission 100" north of the PWN. The Eastern Limb 
is an arced feature that begins at the southern edge and wraps 


6 









The Astrophysical Journal, 789:72 (12pp), 2014 July 1 


Nynka et al. 



Figure 7. Left: NuSTAR 3-6 keV deconvolved image. Intensity profiles were obtained along three lines, shown in yellow, blue, and red, oriented at P.A. = 
0°, 45°, and 90°, respectively. Right: intensity profiles obtained from the deconvolved NuSTAR images. The profiles correspond to the lines of the same color, as 
shown on the left. Scaling was chosen to highlight the Eastern Limb and North Spur. 

(A color version of this figure is available in the online journal.) 


clockwise around to the northern edge, located approximately 
~150" from the center of the plerion. For reference, the outer 
green circle has a radius of 165". Figure 6(b) shows the final 
deconvolved 3-6 keV NuSTAR image, created with 20 iterations. 
The NuSTAR features correspond well to those seen by Chandra , 
allowing us to confidently extend this deconvolution technique 
to higher energies. 

4.2.3. Image Deconvolution: Analysis 

Once we established that the NuSTAR results reproduce 
the Chandra images stably for a range of Lucy-Richardson 
iterations, we applied the same deconvolution to the higher 
energy images. 

Figures 6(b)-(f) show the NuSTAR deconvolved images at 
3-6 keV, 6-10 keV, 10-15 keV, 15-20 keV, and 20-25 keV, 
respectively. One can visually identify the Eastern Limb and 
North Spur at energies as high as 15 keV, as seen in Figure 6(d). 
This is not surprising, since all of the previously reported 
spectral fitting of these two regions have required a non-thermal 
component. This is, however, the first direct measurement 
showing that these features have emission above 10 keV. 

At energies above 15 keV the Limb and Spur become 
very faint. One can still see emission along the eastern edge 
of the Limb in Figure 6(e). We can, however, verify the 
existence of these features above 15 keV by producing intensity 
profiles of the deconvolved images at various azimuthal angles. 
Applying the same one-dimensional (ID) profile method as in 
Section 4.2.1, we obtained profiles along three different radial 
lines, each bisecting G21.5— 0.9 at angles in increments of 45° 
as shown in the left image of Figure 7. The arrows overlaid on 
the lines indicate the direction along which the profiles were 
obtained. 

The intensity profiles themselves for the 3-6 keV image are 
shown on the right in Figure 7. Each profile corresponds to the 
line of matching color overlaid on the deconvolved image on 
the left. One can clearly see a sharp increase corresponding to 
the Eastern Limb in the blue and yellow profiles, while the red 
excess indicates the existence of the North Spur. Similarly, the 
Eastern Limb and North Spur are clearly visible in the profiles 
taken from the images up to energies of 15 keV (see Figure 8). 
The intensity profiles from the 15-20 keV image also show 
statistically significant photon excess from the Limb and Spur, 


confirming the existence of these features to energies as high as 
20 keV. 

5. TIMING SEARCH 

The flux in the 0.2-10 keV band from the wind nebula 
generated by PSR J 1 833 — 1034 completely dominates the pulsar 
itself. The pulsar is barely resolved from the PWN at arcsecond 
resolution by Chandra and completely swamped by the PWN 
emission for other X-ray telescopes. All previous searches for 
X-ray pulsation have been unsuccessful despite extensive X- 
ray data sets collected with sufficient timing resolution. The 
extended energy band of NuSTAR presents a new opportunity to 
further isolate the pulsar signal from the PWN. 

To search for the signal from PSR J1 833 — 1034, we initially 
selected photons in the energy range 10 < E < 70 keV from 
a small source extraction aperture (r < 10") (see Figure 4). 
We searched these photons for significant power around the 
expected period using the two derivative radio timing ephemeris 
presented in Abdo et al. (2010, Epoch 2009) extrapolated 
to the NuSTAR observation epoch. This solution is preferred 
over the slightly updated five derivative model of Ackermann 
et al. (2011) whose extrapolated behavior is not predictive. 
Lacking a coeval ephemeris, it is not possible to maintain phase 
unambiguously over the 206 day gap between the two NuSTAR 
observations. Instead, we separately search data collected in 
2012 and 2013 which span 5.8 and 2.65 days, respectively. 

Taking into account the increased uncertainty in the timing 
parameters for the extrapolated ephemeris, we searched for a 
significant signal over a frequency range of ±3 x ay, where a j 
is the uncertainty in the frequency measurement, oversampled 
by three times the Fourier resolution. We evaluate the power at 
each frequency using the test statistic for n — 1, 2, 3, 5, to 
be sensitive to both broad and narrow pulse profiles, possibly 
single or double peaked. The most significant signal in this 
search range was Z\ — 19.34 and Z\ — 21.05, using 6.4 kcts 
and 13.2 kcts, respectively, for the 2012 and 2013 observations. 
This corresponds to a significance of 0.72 and 0.043 after 
taking into account the number of search trials, 20 and 24, for 
the two observations. We repeated our search for a additional 
combination of energy ranges 10 < E < 20 keV, 20 < E < 
79 keV and aperture size with radius <20" but find no signal with 
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Figure 8. Intensity line profiles and significance of NuSTAR deconvolved images. Excess emission indicating the detection of the Eastern Limb and North Spur is 
visible in all energy bands. The colors correspond to the same angles as indicated in Figure 7. The gray dashed line is the background emission, taken from the east 
side of the PWN, where there is no emission from either the Eastern Limb or the North Spur. 

(A color version of this figure is available in the online journal.) 


a greater significance than the initial 2 a result. We conclude that 
no pulsed X-ray signal in detected from PSR J1833 — 1034 in the 
optimal NuSTAR band and place an upper limit at the 99.73% 
(3a) confidence level on a sinusoidal signal pulse fraction of 
f p ~ 4.2% and 6.1% (including unknown PWN emission) for 
the two observations, respectively. 

NuSTAR is not able to independently measure the spectrum 
of the pulsar. However, Matheson & Safi-Harb (2010) were able 
to isolate and fit the spectrum of a 2" region located at the cite of 
the radio pulsar. Using their best-fit non-thermal spectrum we 
can approximate the contributions of the pulsar and PWN in the 
NuSTAR region. This increases the pulse fraction to f p ~ 19.2% 
and 27.9%, respectively. 

6. DISCUSSION 

6.1. PWN Spectral Break and Softening 

Chandra observations of G2 1.5— 0.9 show spectral softening 
over the PWN, with the photon spectral index ranging from 
r ~ 1.4 at the inner (radius < 5") region to T ~ 2.1 at the 
outer 35"-40" radius annulus. NuSTAR observations confirm 
this by showing spectral softening below 9 keV, as shown in 
Table 1. The broad NuSTAR PSF causes some mixing of the 
spectra in different annuli, and as a result the variation of T with 
radius obtained by NuSTAR is less pronounced than that seen 
by Chandra. NuSTAR finds a spectral break at 9.7 keV in the 
integrated PWN emission. The spectral break is also observed 
in the inner regions with radius <30" and in an annulus from 
35"-40". The break is statistically significant at radii >40". 
Above 10 keV the photon spectral index remains constant. 

There are several possible origins for the spectral break seen 
by NuSTAR. One possibility is that the NuSTAR PSF mixes the 
radially dependent power-law indices seen by Chandra (Slane 
et al. 2000; Safi-Harb et al. 2001), and softens the spectrum 
in such a way as to cause a sharp break. However, simulations 
that fold Chandra maps through the NuSTAR response indicate 
that this is not the case. The effects of the large interstellar 
medium (ISM) extinction (Tsujimoto et al. 2011) again would 
not cause such a sharp, defined energy break. Similarly, the loss 
effect of dust scattering at lower energies is negligible (below 
2%) in the NuSTAR band because of the relatively low column 


absorption. Contributions from the pulsar are also not likely to be 
responsible. Most pulsars have a harder spectrum than the PWN 
they power. While this could in fact cause a sharp transition 
between photon indices, it would cause spectral hardening with 
energy, not the softening seen by NuSTAR. Finally, spectral 
breaks are often attributed to synchrotron cooling, as proposed 
by Tanaka & Takahara (2011). However, the break energy 
of 9 keV would require an unreasonably low magnetic field 
strength of ~6 /xG, compared to the ~300 /xG derived from 
equipartition arguments (Camilo et al. 2006). The prominent 
spectral steepening between radio and X-rays, if interpreted with 
simple cooling models, yields magnetic field strengths ranging 
from 25 /xG (de Jager et al. 2008) to 64 /xG (Tanaka & Takahara 
2011 ). 

It is likely that the spectral break results from physical 
effects, either due to a break in the injected electron energy 
spectrum or due to energy losses due to particle transport in the 
PWN. Pulsars emit pairs of relativistic electrons and positrons, 
which are accelerated at a termination shock near the pulsar 
itself. Downstream of the termination shock the accelerated 
electrons interact with the magnetic field, also produced by the 
pulsar, and subsequently emit synchrotron radiation from the 
radio through gamma energy bands. The injection spectrum can 
therefore shape and influence the spectrum of the synchrotron 
radiation, as noted by Tanaka & Takahara (2011). These authors 
propagate a broken power-law electron spectrum through a time- 
dependent model that includes energy losses due to synchrotron 
radiation, IC scattering, and adiabatic cooling. NuSTAR provides 
confirmation that the relatively simple injection spectrum used 
in Tanaka & Takahara (2011) is not adequate to fit the observed 
X-ray data, as seen in Figure 2. The model, noted as the solid 
black line, has a steep negative slope in the X-ray band and does 
not fit the X-ray spectra obtained from Chandra , XMM-Newton , 
INTEGRAL/ IBIS, and now NuSTAR. 

It is therefore reasonable to explore whether a more complex 
model can explain the 9 keV spectral break. Vorster et al. 
(2013) extend the models with the addition of diffusive losses 
as well as a broken injection spectrum with a discontinuity 
at the break energy. However, both the aforementioned SED 
models do not include the spatial dependence of parameters such 
as the magnetic field of the PWN, which provides additional 
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complexity and can perhaps better fit the X-ray data. We explore 
these effects in the following sections. 


6.2. Physical Conditions Inferred from 
Cooling Scale Length Measurements 

Particle transport in PWNe has long been a matter of dis- 
cussion and study. Unfortunately, few celestial objects are both 
bright enough and large enough to allow distinguishing radially 
dependant features to fit to the various existing models. The Crab 
Nebula, G21.5— 0.9 and 3C 58 are such PWNe. They have been 
frequently observed and analyzed, and have provided valuable 
insights into the physics of magnetohydrodynamic (MHD) flows 
and relativistic shocks that govern the appearance of PWNe 
(e.g., Kennel & Coroniti 1984; Komissarov & Lyubarsky 2004; 
Chevalier 2005; Matheson & Safi-Harb 2010). Extending the 
energy range of analysis will allow for more detailed probing 
of radially dependant physics. This information is important 
both to understand PWNe and to isolate the properties of the 
relativistic shock at the PWN inner edge from the subsequent 
spectral evolution downstream. We shall ask: can the observable 
parameters of G2 1.5— 0.9 from radio to X-ray be reproduced as- 
suming the shock injects only a single structureless power-law 
particle spectrum? 

One such parameter used as a spectral fingerprint of various 
models is the radial dependency of the power-law photon 
index. Numerous observations have confirmed softening of the 
G21.5— 0.9 PWN spectrum with increasing projected radius 
(Safi-Harb et al. 2001; Warwick et al. 2001; Matheson & 
Safi-Harb 2005, 2010; Slane et al. 2000). This is shown at higher 
energies for the NuSTAR observations in Figure 4. The softening 
is associated with synchrotron cooling of the electrons, and the 
corresponding decrease in the maximum emitted X-ray energy 
in the bulk velocity flow downstream of the termination shock. 

An alternative approach is to characterize the variation of 
source size with photon energy using a characterization of the 
“cooling length” L(E ), such as the source FWHM. Since only 
the energy dependence of this length is important for modeling 
the spectral steepening, its precise definition is not important. 
This scale depends only on mapping the total number of photons 
as a function of radius and energy to determine the scale length 
L(E). This is more straightforward than spectral modeling 
because counting statistics are almost always limited, and length 
scale measurements do not require determining a spectrum at 
each radius. We therefore use the cooling length scale, with 
FWHM as its surrogate, as the fundamental diagnostic for 
extracting information about physical conditions. 

We shall also make use of the observed steepening from radio 
to X-rays of G2 1.5— 0.9 through the parameter A = a x — a r . 
Most discussions of PWN physics (e.g., Chevalier 2005) simply 
take this as an intrinsic property, but we shall attempt to explain 
it through evolutionary effects. 

There exist two prominent mechanisms that have been in- 
voked to explain PWNe particle transport: advection and diffu- 
sion. Beginning with Wilson (1972) and Gratton (1972), diffu- 
sion has long been investigated as a cause for the characteristics 
of a PWN. However, the early models of particles propagating 
outward from a central source purely by diffusion, applied ex- 
clusively to the Crab Nebula, have been unable to account for 
the detailed X-ray properties of the Crab, such as its change 
of size with frequency (Ku et al. 1976). Pure advection mod- 
els have also been proposed, such as Rees & Gunn (1974). 


The canonical theory involving pure advection as the method of 
particle transport was presented by Kennel & Coroniti (1984, 
hereafter KC84). When KC84 was used to predict the radial 
behavior of the spectrum, the resultant spectral photon index 
has little to no variation from the center of the nebula outward, 
and begins to vary only toward the PWN periphery (Tang & 
Chevalier 2012, hereafter TC12, Figure 2). This also does not 
match the observed behavior of a slowly steepening X-ray spec- 
trum Matheson & Safi-Harb (2005). 

TCI 2 provide a nuanced approach to particle transport by 
diffusion and present two updated models. They claim that the 
magnetic field is not predominately toroidal far from the termi- 
nation shock, as is often assumed in PWN theory, but has a more 
complex geometry with cross-field transport best described by 
diffusion. Their first model consequently incorporates pure dif- 
fusion only, with both the magnetic field and diffusion coef- 
ficient constant with radius, and synchrotron emission as the 
only loss of energy. TCI 2 argue that such a model better ex- 
plains the radial dependence of T as seen in the Crab, 3C58, and 
G21.5— 0.9, with proper adjustment of the diffusion coefficient. 
Involving complexities such as an energy-dependent diffusion 
coefficient might be more physically reasonable, but the data 
are not sufficiently constraining to distinguish these cases from 
a simple diffusion model. 

While the pure diffusion model of TCI 2 appears to describe 
observations of the Crab and 3C 58 relatively well, we argue 
that it is less appropriate for G2 1.5— 0.9. It does provide a good 
description of T(r) for the Crab and 3C 58 (x 2 ~ 1); however, 
the fit of this model to the G21.5— 0.9 >10 keV data is poor 
(X 2 ~ 3), as seen in Figure 6 in TCI 2. In addition, the ratio of 
advective to diffusive timescales that determines what process 
dominates is radially dependant. Using v a 1/r 2 downstream of 
the PWN termination shock (KC84), and the know relationships 
Ldv oc r/v and Uiff oc r, we obtain /advAdiff ex r. G2 1.5— 0.9 is 
very compact, with a size more than two times smaller than that 
of the other two PWNe, and thus advection is likely to dominate. 

Finally, the advective model provides a good fit to the energy 
dependent cooling length scale (see below). TCI 2 presented 
the radially dependant spectral index T(r) rather than using 
the cooling length scale L(E) as characteristic of their model, 
which makes fitting the model to higher energy NuSTAR data 
too difficult. 

The second model presented by TCI 2 involves a Monte Carlo 
simulation that includes both diffusion and advection transport 
close to and farther away from the pulsar, respectively. This al- 
lows a more complex treatment than previous analytical models 
(Massaro 1985). While likely a more appropriate approach for 
G21.5— 0.9, TC12 only applied this simulation to the Crab and 
3C58. 

KC84 provide a complete advective solution for a steady, 
spherically symmetric wind terminated by a MHD shock. 
Although KC84 is routinely applied to determine quantities such 
as the mean downstream magnetic field, its range of applicability 
is in fact more narrow. KC84 represents an idealized theory, 
suited to the case of constant injection of electrons in a spherical 
geometry, transported outward by pure advection in an ideal 
MHD flow with an ordered, toroidal magnetic field. The model 
does not attempt to reproduce the Crab spectrum from radio to 
X-rays, but only the optical to X-ray portion, with an injection 
spectral index of optically emitting electrons of a 0 = 0.6. The 
predicted value of steepening of a x — a 0 of 0.51 is roughly 
appropriate (so their cooling break is at UV wavelengths). 
This value is fortuitously close to the value of 0.5 for a 
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stationary, homogeneous source. Thus this model also cannot 
reproduce radio-to-X-ray SEDs of many other PWNe, including 
G21.5— 0.9, which have a larger A. As mentioned above, the 
T(r) relationship predicted by the KC84 model is flatter in 
the PWN interior and softens toward the edges, as shown by 
Reynolds (2003) and Tang & Chevalier (2012), while generally a 
gradual variation in T(r) is observed. This has motivated several 
generalizations of KC84. 

The values we obtain from the NuSTAR data for L(E ) and 
a x are clearly inconsistent with the model presented in KC84. 
The gradients implied by the assumptions of KC84, predict 
L(E ) oc E~ l/9 , independent of spectral index, and A = (4+cu)/9 
for the physically important inner flow region. If we attempt 
to apply the KC84 formalism to describe the radio-to-X-ray 
spectrum of G21.5— 0.9 with a r = 0.0 and a x = 0.9, we fail 
on both counts. First, we find m = —0.21 zb 0.01 (Section 4.1). 
Second, we find the observed spectral indices in the radio and 
X-ray bands produce A = a x — a r = 0.9 ± 0.1, instead of 5/9 
from a KC model using the injected spectral index of a = 1. 

Reynolds (2009, hereafter R09) noted that there are a number 
of physical effects not accounted for in KC84 that could produce 
a steeper spectral break than A = 0.5. For example, the magnetic 
field may have a significant radial or turbulent component, or it 
may not satisfy mass conservation due to cloud evaporation, or it 
may not have magnetic flux conservation (e.g., due to magnetic 
reconnection or turbulent amplification). R09 constructed a 
simple model that includes these effects. The model involves 
generating simple scaling relations for the downstream magnetic 
field B , fluid velocity field v , and fluid density p in terms of the 
dimensionless length scale L = r/r 0 , where r 0 is the inner 
injection radius. The non- spherical geometrical effects could 
also be parameterized in terms of jet width w = w 0 L € , where 
6 = 1 corresponds to a conical jet or a section of spherical 
outflow; a confined jet has 6 < 1, while a flaring jet would have 
6 > 1. By assuming B , v and p all vary as power laws in L 
with indices m^, m v , and m p , R09 obtains a series of general 
consistency relations that these indices must satisfy with each 
other and with observable parameters. These variables of m, a, 
and A, the energy scaling of the cooling length, particle injection 
index and spectral index break, respectively, provide constraints 
on the allowed values of 6, m^, m v , and m p . 

The results of R09 can be rewritten in terms of the functional 
form of m, and assuming a r = 0 as observed for G21.5— 0.9: 

A = (— m)( 1 + 26 + m p + m^/e 

independent of m v . Since we observe A = 0.9 and m = —0.21, 
this gives l + m p + = 2.296. The possible solutions 

are thus very restrictive for our measured value of A and 
m. Additionally, if the flow conserves mass (disallowing, for 
instance, mass loading by evaporation of thermal material), 
then the density index m p is linked to the velocity index m v . 
In this case, some solutions are unphysical, such as those 
with m v > 0, corresponding to an accelerating downstream 
fluid. With some judicious rejection of such solutions we 
can draw some interesting general conclusions, based on the 
possible values of the observables, which may guide further 
investigations. 

For a conical or spherical flow, 6 = 1, so either density or 
magnetic-field strength, or both, must rise with radius (since 
either m p or is positive, or both). Mass conservation links 
m p to m v , but if that assumption is abandoned, there is no 
relation and no constraint on the velocity profile — only the 
density profile matters for observable quantities. A steeply 


decelerating flow can produce m p > 0 with mass conservation, 
but mass loading can also do this. Similarly, > 0 can be a 
reasonable outcome of flux nonconservation through processes 
such as reconnection. Our observations require one or both of 
these effects: addition of mass to the flow through some kind 
of evaporation, and increase in magnetic-field strength beyond 
flux freezing. 

A strongly confined jet (6 < 1) can relax some constraints; 
for 6 = 0.3, we only require m p + = —0.32. However, 

even here either mass or flux nonconservation is necessary. One 
can construct constant density solutions ( m p = 0) but such a 
geometry is disfavored due to the high symmetry of G2 1.5— 0.9 
as evidenced in both the broad axial symmetry observed in the 
radio (Furst et al. 1998) and in the soft X-ray band (Safi-Harb 
et al. 2001). 

While Reynolds (2009) expands on the treatment of KC84, 
both fail to reproduce the steadily steepening spectrum with 
radius shown by G21 .5— 0.9 and other PWNe. This shortcoming 
is characteristic of models in which particles are transported 
outwardly purely by advection, so that all particles at a given 
radius have similar ages. To produce steady spectral steepening 
probably requires a mixture of particles of different ages at each 
radius. This could be caused by more complex fluid flow such as 
the back flows found in simulations by Komissarov & Fyubarsky 
(2004), or by particle diffusion. 

We conclude that a model describing both the radio-to- 
X-ray spectrum of G2 1 .5 —0.9 and the size shrinkage with X-ray 
energy we observe can be accommodated in a pure advection 
model requiring the injection of only a straight power-law 
spectrum of electrons, N(E) oc E~ l . However, as with all pure 
advection models, the gradual rather than sudden steepening 
of the spectrum with radius is not reproduced. The viability 
of this explanation for the observed properties of G2 1.5— 0.9 
will depend on whether the addition of diffusion can reproduce 
the gradual steepening while preserving the successes of the 
advection model. 

6.3. The North Spur 

NuSTAR has detected, for the first time, the North Spur 
and Eastern Fimb above 10 keV. Three main theories have 
been proposed to explain the nature of these features: they 
are extensions of the PWN itself, they are limb-brightened 
shock fronts propagating into surrounding ISM and accelerating 
cosmic rays, or they result from an interaction of ejecta with the 
envelope of the progenitor SNR, presumably a Type IIP SN 
(Bocchino et al. 2005). Since the North Spur and Eastern Fimb 
have different spectral and spatial properties, we discuss them 
separately. 

A multi-wavelength analysis is required to get a full under- 
standing of the North Spur. This feature was observed in the 
radio and the soft X-ray, most recently by Bietenholz et al. 
(2011) and Matheson & Safi-Harb (2010), respectively. Bieten- 
holz et al. reported a radio detection of the North Spur, with a 
1.43 GHz flux density of 20.2 ±1.8 mJy, and a FWHM size of 
1 8" x 8". This is notable because it is the only feature, other than 
the PWN itself that is detected in the radio band. The Eastern 
Fimb has no radio emission detected to date. 

Matheson & Safi-Harb (2010) obtained ~580 ks of Chandra 
data, and found that the North Spur has a spectrum comprised of 
non-thermal and thermal components. The thermal component, 
represented by the pshock model, is best fit to temperature 
of kT ~ 0.2 keV and contributes only ~6%-7% to the 
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overall 0.5-8 keV flux. The non-thermal component, however, 
is equally well fit by either the srcut or powerlaw model. 

The model srcut describes the synchrotron emission from 
a homogeneous source consisting of a power-law energy dis- 
tribution of electrons with an exponential cutoff, radiating in a 
uniform magnetic field. The emitted spectrum is a power law 
that steepens slowly above the photon energy corresponding 
to the electron cutoff energy. This slow curvature can mimic 
a steeper power law in a limited energy band. However, in a 
broader energy band the srcut model can fall well below the 
extrapolation of a power law with the same slope at lower en- 
ergies. Since X-ray emission is visible from the North Spur at 
energies as high as 20 keV, the correct spectral model must pro- 
vide a photon flux from 15-20 keV that is statistically higher 
than the background. 

Due to signal-to-noise limitations, we cannot use NuSTAR 
data to spectrally fit the Spur, however, we can use the decon- 
volved images to distinguish between the srcut and powerlaw 
models. We simulated spectra with these models using the 
fakeit command in XSPEC using the model parameters reported 
by Matheson & Safi-Harb (2010). A photon spectral index of 
r = 2.21 was applied to the powerlaw model, while the srcut 
model used a radio index of a = 0.8 and a rolloff frequency of 
Vroiioff = 18 x 10 17 Hz. The NuSTAR response files were based 
on a point source extraction with r = 30". We were thus able 
to obtain fluxes from each respective model from 10-15 keV, 
15-20 keV, and from 20-25 keV. 

Both models have count rates higher than that of the back- 
ground within the 10-15 keV energy band. In the 15-20 keV 
band, the powerlaw model has a count rate seven times higher 
than the background, while the srcut model is only three times 
higher. Finally, within the 15-20 keV band, the count rates are 
4 and 1.5 higher than the background for the powerlaw and 
srcut models, respectively. This implies that the srcut model 
should be marginally visible up to 20 keV, and above 20 keV 
should have a count rate equal to that of the background. This 
matches well with what is seen in the NuSTAR images. The 
powerlaw model, however, should be detectable at energies as 
high as 25 keV. 

If the North Spur were an extension of the PWN, it would have 
a spectral photon index similar to that of the PWN itself. While 
this is true, the analysis above indicates that the North Spur is not 
described by a powerlaw model extending to higher energies. If 
it were, the North Spur would easily be detectable at energies as 
high as 25 keV. However, the NuSTAR images with a combined 
exposure of 281 ks do not detect any statistical emission above 
20 keV, indicating that srcut is the more plausible spectral 
model for the North Spur. 

It is possible that the North Spur results from the interaction 
of the inner SN ejecta with the H-envelope of the progenitor 
(Bocchino et al. 2005; Matheson & Safi-Harb 2010). This 
is supported by the thermal component in the spectral fit of 
Matheson & Safi-Harb (2010), a pshock model with solar 
abundances, low temperature of kT — 0.21 ± 0.4 keV, and low 
ionization timescales. This is also supported by the morphology 
of the North Spur itself. With the inclusion of projection effects, 
the North Spur is located between —75" and the edge of the SN 
shell at -120". 

6.4. The Eastern Limb: The Shell of G2 1.5— 0.9 

X-ray emission from the shell of G21.5— 0.9 is clearly 
visible in the NuSTAR image. Our results confirm the existence 


of this shell, revealed by Chandra (Matheson & Safi-Harb 
2010) in both their image and in an extracted shell spectrum 
(up to —6 keV), but also hinted at in earlier XMM-Newton 
data (Bocchino et al. 2005). The XMM-Newton observations 
reveal evidence of the shell in a —2-8 keV energy band 
image, after careful subtraction of a modeled dust scattering 
component below —5 keV. The NuSTAR detection extends up 
to much a higher energy of —20 keV. The morphology of the 
NuSTAR emission is striking in its similarity to the Chandra and 
especially the XMM-Newton image. Emission is detected from 
position angle (P.A.) -180° to -300° in the 6-10 keV image 
(P.A. = 0 at north, positive clockwise), with the extent shrinking 
as the energy increases until it is visible mainly in the east and 
north at the highest energies. This is consistent with the intensity 
distribution with position angle seen in the lower energy image. 

A question unresolved by previous X-ray observations is 
whether the shell emission is thermal or non-thermal. The 
extended energy response of NuSTAR can be exploited to answer 
this question. Matheson & Safi-Harb (2010) found that the 
spectrum of the Eastern Limb could be characterized equally by 
four distinct models. A thermal fit to the pshock model provided 
a temperature of kT — 7.5 keV while a non-thermal powerlaw 
model obtained a spectral photon index of T = 2.13. Two srcut 
models were also well fit to the Eastern Limb spectrum. Ideally 
srcut requires both a radio flux density and radio spectral 
index for the shell, but the Eastern Limb has not been detected 
in the radio. Therefore, two values of the radio index a that 
are reasonable for a SN shell were chosen (< a = 0.3/0. 5) while 
all other parameters were left free to vary. Care was taken to 
ensure the best-fit radio flux was below the upper limit reported 
by Bietenholz et al. (201 1). 

As with the North Spur, we extrapolated the four Chandra 
spectra for the Eastern Limb into the NuSTAR band. We created 
an effective area file for an extended source shaped like the 
Eastern Limb, then used the fakeit command in XSPEC to 
simulate spectral data. The thermal pshock model predicted 
a shell flux which would not produce the X-ray emission 
seen by NuSTAR at —15-20 keV. In contrast, the three non- 
thermal models produced X-ray fluxes consistent with imaging 
of the Eastern Limb by NuSTAR , although the srcut and 
powerlaw models could not be distinguished from each other. 
Nonetheless, the NuSTAR observations firmly establish the non- 
thermal nature of the shell X-ray emission. 

The detection of a non-thermal shell up to quite high 
X-ray energy in G21.5— 0.9 is interesting, and is in contrast 
to observations of other Crab-like SNRs. Recently shells have 
been detected in 3C58 (Gotthelf et al. 2007) and G54. 1+0.3 
(Bocchino et al. 2010), and a clear, detached shell of emission 
in Kes 75 (Helfand et al. 2003). However, the 3C58 shell is 
clearly thermal, with no sign of a non-thermal component. The 
shell of G54. 1+0.3 can be fit with both thermal and non-thermal 
models; however, a thermally emitting shell seems much more 
likely. Assuming thermal emission, Bocchino et al. (2010) were 
able to use PWN-SNR evolutionary models to obtain an age 
consistent with the characteristic age derived from pulsar obser- 
vations, obtain the proper SN-PWN radius ratio, and predict that 
the reverse shock has not encountered the PWN yet, consistent 
with other observations. The ages of 3C58 and G54. 1+0.3 are 
—3000-5000 yr and —1800-3300 yr, respectively. G21.5— 0.9 
is thus unique in that it is much younger (—290-1000 yr), and 
potentially has higher forward shock speed, both of which could 
lead to the detectable non-thermal shell. The Crab Nebula itself, 
of comparable age to G2 1.5— 0.9 but much closer, still shows 
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no non-thermal shell, presumably due to a very low interstellar 
medium density. The detection of this G21.5— 0.9 non-thermal 
shell at quite high X-ray energy is thus likely due to its younger 
age compared to these other Crab-like SNRs. 

7. SUMMARY 

We have presented an analysis of a 281 ks NuSTAR observa- 
tion of G2 1.5— 0.9 to probe the spatial and spectral characteris- 
tics revealed by high-energy X-ray emission. NuSTAR ' s broad 
energy band makes it uniquely suited to analyze not only the 
PWN but observe the characteristics of the SNR shell as well. 
The spectrum taken from the entire remnant is described by a 
best- fit broken power law with a break energy of ~9 keV. This is 
the first instance where a single instrument was able to capture 
this break. The addition of the NuSTAR spectrum to SED models 
from Tanaka & Takahara (201 1) produces a poor fit to the X-ray 
data. This suggests that further modeling is required: more com- 
plex electron injection spectra, additional loss processes such as 
diffusion, or radial dependence of the PWN parameters. 

Spectra extracted from various radial annuli were also fit with 
both an absorbed power-law and an absorbed broken power- 
law model. The two regions r ^ 30" and r = 30"-60" are 
statistically better fit with a broken power law with £break ~ 

9 keV, while the regions with radii r > 60" are best fit with a 
single power law. We observe spectral softening of the spectral 
index below the spectral break, while the spectral index above 
the break is constant within uncertainties. 

Image analysis allows us to measure the energy-dependent 
cooling length scale and fit the relationship with a power- 
law model of L(E) ex E m . This yields an index of m = 
—0.21 zb 0.01. Incorporating this with the spectral indices in 
both the radio and X-ray bands, we are able to systematically 
approach the equations of Reynolds (2009) and inspect the 
parameter space for physically consistent solutions. We found 
that, for a conical jet or spherical outflow, the most reasonable 
solutions do not conserve magnetic flux but do conserve mass, 
indicating turbulent magnetic field amplification. The bulk 
velocity decelerates steeper than that predicted by KC84. 

We detect the Eastern Limb and North Spur at energies above 

10 keV for the first time. A deconvolution method provides 
clear evidence of emission from the North Spur up to 20 keV. 
Extrapolation of the spectral fits obtained by Chandra show 
that the srcut model is favored over the powerlaw model. 
This further solidifies the assumption that the North Spur is an 
interaction of the SN ejecta with the remnant. We also detect the 
Eastern Limb up to 20 keV. We have confirmed the existence 
of non-thermal emission from the Limb, and conclude this 
faint feature is the SN shell of G21.5— 0.9. We are unable to 
distinguish between the non-thermal models fit to the Eastern 
Limb by Chandra. 
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